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High resolution total and indium differential atomic pair distribution functions (PDFs) for 
lno.5Gao.5As alloys have been obtained by high energy and anomalous x-ray diffraction experi- 
ments, respectively. The first peak in the total PDF is resolved as a doublet due to the presence of 
two distinct bond lengths, In-As and Ga-As. The In differential PDF, which involves only atomic 
pairs containing In, yields chemical specific information and helps ease the structure data interpre- 
tation. Both PDFs have been fit with structure models and the way in that the underlying cubic 
zinc-blende lattice of Ino.5Gao.5As semiconductor alloy distorts locally to accommodate the distinct 
In-As and Ga-As bond lengths present has been quantified. 
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I. INTRODUCTION 

Ternary semiconductor alloys, such as In^jGai-xAs, 
are technologically important because they allow the 
semiconductor band-gap to be varied continuously be- 
tween the band-gap values of the end members, GaAs 
and InAs, by varying the composition, x |^,|^. This has 
made the technological characteristics, physical proper- 
ties and structure of Ina;Gai_a;As alloys a subject of nu- 
merous experimental and theoretical investigations. It 
has been found that the lattice parameters of the alloys 
well interpolate between those of the end members which 
is consistent with the so-called Vegard's law. According 
to Vegard's law, the structure of alloys adjusts itself so 
that the individual bond lengths are taking equal, com- 
positionally averaged, values and the bond angles remain 
unperturbed from their ideal values for any alloy compo- 
sition. For this reason, all structure dependent prop- 
erties of Ina;Gai_a;As alloys, including electronic band- 
structure, are often calculated within the virtual crystal 
approximation (VGA) . In this approximation the alloy is 
assumed to be a perfect crystal with all atoms sitting on 
ideal lattice sites and site occupancies given by the av- 
erage alloy composition. Both GaAs and InAs have the 
zinc-blende structure (F43m) where In and Ga atoms 
occupy two interpenetrating face-centered-cubic (fee) lat- 
tices and are tetrahedrally coordinated to each other . 
Accordingly, the VGA assumes that Ina;Gai_2;As alloys 
have the same zinc-blende structure and, furthermore, 
the first neighbor interatomic distances (i.e. In-As and 
Ga-As bonds lengths), bond ionicity, atomic potential 
etc. take some average values for any composition x. 
However, both extended x-ray absorption fine structure 
(XAFS) experiments Q and theory Q have shown that 
Ga-As and In-As bonds do not take some average value 



but remain close to their natural lengths h°Qa-As — 
2.437 A and L°jn-As = 2.610 A in the alloy. This behav- 
ior is close to the so-called Pauling model Q which as- 
sumes that the bond length between a given atomic pair 
in an alloy is more or less a constant, independent on 
composition x. This finding shows that the zinc-blende 
lattice of Ina;Gai_a;As is significantly deformed to accom- 
modate the two-distinct Ga-As and In-As bond lengths 
present. Also, the deformation seems to be confined lo- 
cally since the average crystal symmetry and structure is 
still of the cubic zinc-blende type as manifested by the 
Bragg diffraction patterns. 

It is well recognized that local structural distortions 
and associated fiuctuations in atomic potentials signifi- 
cantly affect the properties of materials and, therefore, 
should be accounted for in more realistic theoretical cal- 
culations 0. Thus it seems there is a clear need for 
more detailed determination of the local structure of 
InxGai_a;As semiconductor alloys which then may be 
used as an improved quality structure input to theoret- 
ical calculations. A number of authors P,p|J^-[l0|] have 
already proposed model structures for these alloys but 
there has been limited experimental evidence to date. 
This prompted us to undertake an extensive experimen- 
tal study of the local atomic structure of In-Ga-As semi- 
conductor alloys. 

The technique of choice for studying the local structure 
of semiconductor alloys has been XAFS 11-1^]. How- 
ever, XAFS provides information about the immediate 
atomic ordering (first and sometimes second coordina- 
tion shells) and all longer-ranged structural features re- 
main hidden. To remedy this shortcoming we have taken 
the alternative approach of obtaining atomic pair distri- 
bution functions from x-ray diffraction data. 

The atomic pair distribution function (PDF) is the in- 
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stantaneous atomic number density-density correlation 
function which describes the atomic arrangement in ma- 
terials [Q. It is the sine Fourier transform of the ex- 
perimentally observable total structure factor obtained 
from a powder diffraction experiment. Since the total 
structure factor, as defined in Ref. 14, includes both 
the Bragg scattered intensities and the diffuse scattering 
part of the diffraction spectrum its Fourier associate, the 
PDF, yields both the local and average atomic structure 
of materials. By contrast an analysis of the Bragg scat- 
tered intensities alone, by a Rietveld type analysis for 
instance, yields the average crystal structure only. De- 
termining the PDF has been the approach of choice for 
characterizing glasses, liquids and amorphous materials 
for a long time [ll7| . However, its wide spread applica- 
tion to crystalline materials, where some deviation from 
the average structure is expected to take place, has been 
relatively recent . The present study is a further step 
along this line. 

Very high real space resolution is required to dif- 
ferentiate the distinct Ga-As and In-As bond lengths 
present in In^-Gai-^As alloys. High resolution is at- 
tained by obtaining the total structure factor S(Q), where 
Q = 47rsin0/A is the magnitude of the wave vector, to 
very a high value of Q ( Q > 40 A^^). Here, 29 is the an- 
gle between the directions of the incoming and outgoing 
radiation beams and A is the wavelength of the radiation 
used. Recently, we carried out a high energy (60 keV; 
A = 0.206 A) x-ray diffraction experiment and succeeded 
in obtaining PDFs for Inj;Gai_2;As crystalline materials 
(x = 0, 0.13, 0.33, 0.5, 0.83, 1) of resolution high enough 
to differentiate Ga-As and In-As first atomic neighbor 
distances present An analysis of the experimen- 

tal data (see Fig. 4 in Ref. |l^) showed that the local 
disorder in Ina;Gai_a;As materials peaks at a composi- 
tion X = 0.5. This observation suggested Ino.sGao.s as 
the most appropriate candidate for studying the effect of 
bond-length mismatch on the local structure of In-Ga- 
As family of semiconductor alloys. An important detail 
of the high energy experiments carried out is that low 
temperature (10 K) was used to minimize the thermal 
vibration in the samples, and hence increase the sensi- 
tivity to intrinsic atomic displacements. This left open 
the question about the impact of temperature on the lo- 
cal structure of luxGai-xAs alloys and necessitated the 
carrying out of an complimentary experiment at temper- 
atures considerably higher than 10 K. To partially com- 
pensate for the inevitable loss of resolution from the ther- 
mal broadening of atomic pair we carried out an anoma- 
lous scattering experiment and determined the In differ- 
ential PDF at room temperature. 

In the present paper we report the high energy low 
temperature and anomalous scattering (In edge) room 
temperature experiments on the Ino.5Gao.5As. The ex- 
perimental total and In differential PDFs have been fit 
with structure models and the way in which the zinc- 
blende lattice locally distorts to accommodate the two 
distinct Ga-As and In-As bonds present has been quan- 
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II. EXPERIMENTAL DETAILS 

A. Sample preparation 

The Ino.5Gao.5As alloy was prepared by mixing reagent 
grade GaAs and InAs powders in the proper amounts. 
These were sealed under vacuum in quartz tubes. The 
powders were heated above the liquidus and held for 3 
hours to melt them, followed by quenching into cold wa- 
ter. The resulting inhomogeneous alloys were ground, 
resealed in quartz tubes under vacuum, and annealed at 
a temperature just below the solidus of the alloy for 72- 
96 hours. This procedure was repeated until the samples 
were homogeneous as determined from an x-ray diffrac- 
tion measurement. The sample for high-energy x-ray 
diffraction measurements was a thin layer of the powder 
held between Kapton foils. The thickness of the layer 
was optimized to achieve a sample absorption ^t ~ 1 for 
the 60 KeV x-rays. A standard sample holder with a 
cavity of rectangular shape (2 cm x 4 cm) and depth of 
0.5 mm was used with the anomalous x-ray diffraction 
experiments. The powder was loaded into the cavity to 
avoid any texture formation and its extended surface left 
openly exposed to the x-ray beam. 



B. High-energy x-ray diffraction experiments 

The high resolution total-PDF measurements and data 
analysis has been reported elsewhere |19 2^. Here, the 
experiment procedures and data analyses employed are 
reported in some more detail. The experiments were car- 
ried out at the A2 24 pole wiggler beam line at Cornell 
High Energy Synchrotron Source (CHESS). All measure- 
ments were done in a symmetrical transmission geometry 
at 10 K. The polychromatic incident beam was dispersed 
using a double crystal Si(lll) monochromator and x-rays 
of energy 60 keV (A = 0.206 A) were employed. An in- 
trinsic Ge detector coupled to a multi-channel (MCA) 
analyzer was used to detect the scattered radiation. By 
setting proper energy windows we were able to extract 
the coherent component of the scattered x-ray intensi- 
ties during data collection. The diffraction data were 
collected in scanning at constant AQ steps of 0.02 A"-'^. 
Several runs were conducted and the resulting spectra 
averaged to improve the statistical accuracy and reduce 
any systematic error due to instability in the experimen- 
tal set-up. The diffraction data were smoothed using the 
Savitzky, Golay procedure [ pl| . The procedure was tuned 
in such a way that each data point gained or lost only 
one Poisson counting standard deviation in the smooth- 
ing process. The data were normalized for flux, and cor- 
rected for background scattering and experimental effects 
such as detector deadtime and absorption. The part of 
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FIG. 1. Total (lower part) and In differential (upper part, 
offset for clarity) reduced structure factors for luo.sGao.s. 



Compton scattering at low values of Q not eliminated 
by the preset energy window was removed analytically 
applying a procedure suggested by Ruland |22|. The re- 
sultant intensities were divided by the square of the aver- 
age atomic form factor for the sample to obtain the total 
structure factor S{Q), 



S{Q) = 1 



(1) 



where F°'' is the coherent part of the total diffraction 
spectrum; Ci and fi{Q) are the atomic concentration and 
scattering factor of the atomic species of type i [i = 
In,Ga,As), respectively llSyiTf] . All data processing pro- 
cedures were done with the help of the program RAD 
||. The reduced structure factor F{Q) = Q[S{Q) - 1] 
is shown in Fig. |l|. As can be seen in the figure the F(Q) 
data are terminated at Qmaa; = 45 beyond which, 
despite the high intensity synchrotron source employed 
and extra experimental data averaging applied, the sig- 
nal to noise ratio became unfavorable. It should be noted, 
however, that this is a very high value of the wavevec- 
tor for an x-ray diffraction measurement; for comparison 
Qmaa; achicvcd with a conventional source such a Cu an- 
ode tube is less than 8 A^^. The corresponding reduced 
atomic distribution function, G{r), obtained through a 
Fourier transform 

G(r) = ATrr[p{r) - po] 




FIG. 2. The reduced total (full line) and In differential 
(symbols) atomic pair distribution functions for Ino.5Gao.5As. 
The In-differential PDF is significantly broader than the to- 
tal-PDF because Qmax (as shown in Fig. 0), and therefore the 
real-space resolution of the measured PDF, is much lower. 



= (2/^) 



Q=0 



Q[S{Q)-1] sin Qr dQ, 



(2) 



is shown in Fig. 2. Here p{r) and po are the local and 
average atomic number densities, respectively, and r the 
radial distance. It should be noted that no modification 
function, i.e. additional damping of the S(Q) data at 
high values of Q, was carried out prior to Fourier trans- 
formation. This ensures that the data have the highest 
resolution possible but it results in some spurious ring- 
ing in G(r). If Qmax is high enough the ringing is small 
(on the level of the random noise), and in any case it is 
properly modeled by convoluting G(r) with a Sine func- 
tion |p3] which we do in all our modeling. 



C. Anomalous x-ray diffraction experiments at the 
In Edge 

It is well known that a single diffraction experiment 
on an rt-component system yields a total structure factor 
which is a weighted average of n(n+ l)/2 distinct partial 
structure factors, i.e. 



5(Q)=^my(Q)5,,(Q), 



(3) 
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where Wij{Q) is a weighting factor and Sij{Q) the par- 
tial structure factor for the atomic pair {i, j), respectively 
psf . The corresponding total PDF, too, is a weighted av- 
erage of n(n+ 1)/2 partial pair correlation functions. For 
a multi component system like Ino.5Gao.5As it is there- 
fore difficult to extract information about a particular 
atomic pair from a single experiment. The combination 
of a few conventional experiments (let say a combination 
of neutron, x-ray and electron diffraction experiments) 
or the application of anomalous x-ray scattering allows 
the determination of chemical-specific atomic pair distri- 
butions. 

We briefly describe the use of anomalous scattering to 
obtain chemical specific PDFs . If the incident x-ray 
photon energy is close to the energy of an absorption edge 
of a specific atom in the material, the atomic scattering 
factor should be considered a complex quantity depen- 
dent on both wavevector Q and energy E 



fiQ,E) = UQ) + fiE)+ir{E), 



(4) 



where fo{Q) is the usual atomic scattering factor and 
/' and /" are the anomalous scattering terms depend- 
ing on the x-ray photon energy E. The imaginary term, 
/", is directly related to the photoelectric absorption co- 
efficient and it is small and slowly varying for E below 
the edge, rises sharply at the edge, and then gradually 
falls off. /' has a sharp negative peak at the edge with a 
width which is typically 40-80 eV at half maximum and 
is small elsewhere . This behavior can be clearly seen 
in Fig. 1^. The anomalous scattering technique takes ad- 
vantage of the fact that /' and /" for a particular atomic 
species change rapidly only within ~ 100 eV of the re- 
spective absorption edge and that the characteristic ab- 
sorption edges for different atomic species are separated 
by several hundreds of eV. Then a difference of two sets 
of diffraction data taken at two slightly different ener- 
gies below an absorption edge of a particular element 
will contain only a contribution of atomic pairs involving 
that element. Accordingly, one can define a differential 
structure factor, DSF, as follows: 



DSF = 



IjE,) - I{E2) - - ^c,p{E2)\ 



(5) 



where Ei and E2 are the two photon energies used P6[ . 
The differential PDF, which gives information about 
the atomic distribution around the anomalous scatter- 
ing atoms, is calculated analogous to Eq. (3) with 
S(Q) replaced by the DSF. Several experiments have al- 
ready demonstrated the usefulness of anomalous scatter- 
ing techniques in studying the local atomic ordering in 
both disordered and crystalline materials p7|-p9[ . A pre- 
cise knowledge of the anomalous scattering terms is, how- 
ever, a prerequisite for the interpretation of anomalous 
scattering experiments. Unfortunately, theoretical mod- 
els are not capable of providing precise enough values for 
/' and /" in the vicinity of absorption edges. That is 
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FIG. 3. Energy dependence of the real, /', and imagi- 
nary, /", anomalous scattering terms of atomic x-ray scat- 
tering factor of In. (Theoretical data - symbols. The dotted 
line through the symbols is a guide to the eye. Experimental 
data - full line). The two energies below the In edge used in 
the present anomalous scattering experiment are marked by 



why anomalous scattering experiments usually involve a 
complimentary determination of /' and /". It is most 
frequently done by measuring the energy dependence of 
/" in the vicinity of the absorption edge and a subsequent 
determination of /' through the so- called dispersion re- 
lation iSOll: 



nE,) = (2/^) 



[El^E^- 



dE 



(6) 



The same strategy was adopted in the present anomalous 
diffraction experiments. These were carried out at the In 
edge which is the highest energy (~ 27.929 KeV) edge ac- 
cessible in In-Ga-As system. The experiments were car- 
ried out at X7A beam line at the National Synchrotron 
Light Source at Brookhaven National Laboratory. Two 
energies, one just below (27.889 KeV) , and the other 
few hundred eV (27.629 KeV) below the In edge were 
used. The raw data are shown in Fig. ^. A Si channel- 
cut monochromator was used to disperse the white beam. 
The monochromator was calibrated by measuring the ab- 
sorption edge of indium of high purity. The scattered 
x-rays were detected by a Ge solid state detector cou- 
pled to a multi-channel analyzer. Few energy windows, 
covering several neighboring channels, were set up to ob- 
tain counts integrated over specific energy ranges during 
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FIG. 4. Coherently scattered intensities of an 

Ino.5Gao.5As sample measured at two different energies just 
below the In edge (dots and solid line, respectively). The dif- 
ference, A7, (offset for clarity) between the two data sets is 
given in the lower part. 



the data-collection. These energy windows covered the 
coherent intensity; coherent, incoherent and In K[3 fluo- 
rescence intensities all together; In Ka fluorescence; As 
Ka fluorescence, and a window covering the entire energy 
range which integrates the total scattered intensity in the 
detector. Integrated counts of these ranges were collected 
several times and then averaged to improve the statisti- 
cal accuracy. The data were corrected for detector dead 
time, Compton and background scattering, attenuation 
in the sample and residual In KP fluorescence which is 
not possible to be well separated from the coherent com- 
ponent of scattered intensities. In K(3 was determined by 
monitoring the In Ka signal and multiplying it by the 
ratio of In K(3 to In Ka output, which was experimen- 
tally determined by a complimentary experiment carried 
out well above the In edge (~ 29 keV). 

By taking the difference between the two data sets, 
as shown in Fig. ^, all terms that do not involve In were 
eliminated, because only In scattering factor changed ap- 
preciably in the energy region explored while the scatter- 
ing factors of Ga and As remained virtually the same, and 
the In DSF was obtained. 

The unknown anomalous scattering terms of In, in- 
volved in Eg. 5, were determined in the following way: 
In fluorescence yield was detected by scanning over a wide 
range across the In edge. The curve was matched to the 



theoretical estimates of Chantler for /" and so the 
fluorescent yield converted to /" data. /' was calculated 
from these /" data via the dispersion relation as given 
in Eq. 6. The anomalous scattering terms of In, resulted 
from the present experiments, are given in Fig. ^. As 
one can see in the figure, in the vicinity of In edge /' and 
/" change sharper than theory predicts. We determined 
the following values for /' and /" for the two energies 
employed: /' = -3.89 and /" = 0.637 at £: = 27.629 
keV; /' = -6.148 and /" = 0.826 at E=27,889 keV. The 
use of the experimentally determined but not the theo- 
retically predicted values of /' and /" turned out to be 
rather important in obtaining differential structure data 
of good quality. The In-DSF and differential PDF for 
Ino.5Gao.5As alloy are shown in Figs. |l| and 2, respec- 
tively. The In-DSF appears broader primarily because 
Qmaxi and therefore the resolution of the measurement, 
is much lower than is the case for the total-PDF mea- 
surement (as is obvious in Fig. |l|). There is an additional 
broadening of this peak because the In-DSP data were 
collected a room temperature instead of lOK, but this is 
expected to be small. 



III. RESULTS 

As can be seen in Fig. |l| significant Bragg scatter- 
ing (well defined peaks) is present up to approximately 
15 A"-'^ in the In difference and total structure factors 
of Ino.5Gao.5As alloy. At higher wavevectors only an os- 
cillating diffuse scattering is evident. This implies that 
although the sample still has a periodic structure it con- 
tains considerable local displacive disorder. The disorder 
is due to the mismatch of Ga-As and In- As bond lengths 
clearly seen as a split first peak in the total PDF of Fig. 2 
Jigf . Also shown in Fig. 2 is the In difference PDF which 
has a single first peak well lining up with the higher-r 
component of the first peak in the total PDF. Since the 
In differential PDF contains only atomic pairs involving 
In its first peak can be unambiguously attributed to In- As 
atomic pairs. This allows us to identify the two compo- 
nents of the first peak in the total PDF as being due to 
Ga-As and In- As atomic pairs, respectively. According to 
the present high resolution x-ray diffraction experiments 
Ga-As and In-As bond lengths in the Ino.5Gao.5As al- 
loy are 2.455(5) A and 2.595(5) A at 10 K, respectively. 
In the present anomalous diffraction experiments In-As 
bond length is 2.615(5) A at room temperature. The ob- 
served elongation of the In-As bond with temperature is 
due to the usual thermal expansion observed in materi- 
als. We note that the present PDF-based results are in 
rather good agreement with the XAFS results of Mikkel- 
son and Boyce |4j for Ga-As and In-As bond lengths in 
Ino.5Gao.5As. 

An inspection of the experimental PDFs in Fig. 2 (see 
also Figs. H and ^ shows that the nearest atomic neigh- 
bor peak is the only one which is relatively sharp. Start- 
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ing from the second-neighbor peak onwards all atomic- 
pair distributions (PDF peaks) show significant broad- 
ening. The observation shows that the bond-length mis- 
match gives rise to a considerable deformation of the 
underlying zinc-blende lattice of Ino.5Gao.5As alloy. To 
quantify this deformation we explored a few structure 
models as follows: 



A. Supercell model based on the Kirkwood potential 

It was previously shown that a 512 atom supercell for 
the alloy structure, based on the zinc-blende unit cell but 
with In and Ga randomly arranged on the metal sublat- 
tice and atomic positions relaxed using the Kirkwood po- 
tential 1^ , explains well the high-resolution total-PDFs 
In addition to this high spatial resolution PDF we 
now have a PDF which is chemically resolved. We are 
interested to know whether this supercell model is still 
sufficient for describing these new data. The model has 
been described in detail elsewhere |]3^ . 

In the present modeling we used the same force con- 
stants aij and Pij that were selected to fit the end mem- 
bers GaAs and InAs jl^,^^. Using the relaxed atomic 
configuration a dynamical matrix has been constructed 
and the eigenvalues and eigenvectors found numerically. 
From this the Debye- Waller factors for all the individual 
atoms in the supercell have been determined. The PDF 
of the model was then calculated using a Gaussian broad- 
ening of the atomic-pair correlations to account for the 
purely thermal and zero-point motion. The width of the 
Gaussians was determined from the theoretical Debye- 
Waller factors [Q. In addition, the calculated PDF was 
convoluted with a Sine function to account for the trun- 
cation of the data at Q„ 



l33[ . A comparison between 
the model and experimental results is shown in Fig. |^. 
The agreement with both the high spatial resolution data 
(Fig. 1(a)) and with the differential PDF (Fig. |(b))is 
clearly very good. It has already been demonstrated 
that the Kirkwood-based model well reproduces the dis- 
placements of As and metal atoms in In-Ga-As alloys 
extracted from a model independent analysis the PDF 
peak widths. Thus the present and previous results sug- 
gest that the Kirkwood-based model is a good starting 
point for any further calculations requiring good knowl- 
edge of the local structure of the Ino.5Gao.5As alloy. 



B. Refinement of chemically resolved differential and 
Spatially resolved Total PDF 

In this paper we present, for the first time, both high- 
resolution total- and chemically resolved In-partial PDFs 
for Ino.5Gao.5As. In the previous Section we showed that 
the data are consistent with a supercell model based on 
the Kirkwood potential. However, in addition we would 
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FIG. 5. Experimental (open circles) and model (solid line) 
PDFs for Ino.5Gao.5As alloy, (a) total PDF; (b) In differen- 
tial PDF. Model PDFs are based on Kirkwood-type potential 
minimization. The residual difference between the model and 
experimental data is given in the lower part. 



like to extract structural information from the data with- 
out recourse to potential-based models which can be used 
to compare with other experimental results |^,^ [ 
and theoretical predictions PJ^-pO[|. 
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To do this we have constructed the simplest possible 
model that was still consistent with the data, and we have 
refined it using the PDF profile-fitting program PDF- 
FIT . We have fit to both the high-spatial resolution 
total-PDF and the chemically resolved differential PDF 
data at the same time which resulted in equivalent atomic 
displacement parameters being refined. 

The model is based on the 8-atom cubic unit cell of 
the zinc-blende structure. The split nearest-neighbor 
peak in the total PDF, and the shifted nearest-neighbor 
peak in the In-differential PDF, both require that defi- 
nite static displacements of fixed length be incorporated 
in the model. In addition, the shifted position of the 
nearest neighbor peak in the In-differential PDF requires 
that a model be constructed which has a definite chemi- 
cal species on specific sites, i.e., goes beyond the virtual 
crystal approximation. In simple 8-atom cubic unit cell 
this de facto leads to a chemically ordered model that 
is not observed in the real alloy and which, furthermore, 
does not sample all of the possible chemical environments 
for As in the random alloy HJl^,^. Nonetheless, this is 
the minimal model which can be successfully refined to 
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the experimental data to extract information about local 
atomic displacement amplitudes. 

In this model the four metal sites are populated with 
two In and two Ga ions. Static displacements of As and 
metal ions were then allowed. The model was constrained 
so that all four As sites had the same displacement am- 
plitude. The metal sites were likewise constrained to be 
displaced by the same amount as each other, but the 
metal site displacement was independent of the As sub- 
lattice displacement. The directions of the displacements 
were also constrained to be along (111) type directions. 
The choice of which of the 8 possible (111) directions 
was determined by the chemical environment. A model 
with (100) type displacements was less successful at rec- 
onciling the sharp first peak and broad later peaks in 
the PDF. This is discussed in more detail later. We 
call these the "discrete" displacements. In addition to 
the discrete atomic displacements, atomic-displacement- 
parameters (thermal factors) and lattice parameters were 
refined. The thermal factors contain both static and dy- 
namic disorder. These atom displacement distributions 
we refer to as "continuous" to differentiate from the dis- 
crete displacements described above. Finally, the near- 
est neighbor peak was sharpened with respect to the rest 
of the PDF using a sharpening factor. This accounts 
for the highly correlated nature of the displacements of 
near-neighbor atoms. The resulting fit to the data is 
shown in Fig. |^. 

The values we refine are as follows: the discrete 
displacements on the As and Metal sublattices are 
0.133(5) A and 0.033(8) A respectively. These val- 
ues are independent of temperature. The continuous- 
displacement amplitudes are {u\^ ) = 0.00814(12) A2 and 
0.00373 A^ for the As and metal sublattices 



0.0135(15) A2 and 
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at T = 10 K and {u\^ 

0.010(2) A^ , respectively, at room temperature. These 
compare with literature values of {u\^) = 0.0015(8) A^ 
and {u1[) = 0.0017(9) A^ for the end-member com- 
pounds at 10 K |l9|, and {u\^) = 0.00716(5) A^ and 
(""m) = 0.009 A^ at room temperature for As and 
the metal site, respectively. 

The discrete displacements obtained from the fits are 
illustrated schematically, by a projection of a fragment 
of the structure down the [010] direction, in Fig. |7| This 
shows that both the discrete and continuous displace- 
ments (indicated schematically by the size of the circles 
representing the atoms) are larger on the As than the 
metal sublattice. The size of the circles representing the 
continuous-displacements have been exaggerated. 



IV. DISCUSSION 

Existing experimental results that characterize the 
structure of semiconductor alloys beyond the average 
structure include XAFS_|0,§,|5| , ion-channeling Q, x- 
ray diffuse scattering [|36[, and Raman scattering [ p7[ . 



O 




FIG. 6. Experimental (open circles) and model (solid 
line) PDFs for Ino.5Gao.5As alloy. (a) total PDF; (b) 
In differential PDF. Model PDFs are based on a distorted 
zinc-blende-type structure. The residual difference between 
the model and experimental data is given in the lower part. 



The XAFS results clearly show that short and long near- 
neighbor bonds exist which are from Ga-As and In-As 
neighbors respectively. The bond-length distribution is 
not recovered with great accuracy and there is limited 
information available on higher neighbor pairs; however, 
the indication is that the atom-pair separations return 
quickly to the virtual crystal values with increasing pair- 
separation, r. This implies significant distortions to 
the crystal structure. Indeed, a correct analysis of the 
phonons in Raman spectra from Gai_a;In2;As required 
significant structural distortions pTf . Our current and 
earlier |l^ PDF results bear out all these observations. 

The discrete displacements refined in our (111) dis- 
placed model are primarily determined by the splitting, 
and displacement, of the first peak in the total and In- 
differential PDF's, respectively: this sharp feature in the 
PDF is very sensitive to the amplitude of the discrete 
displacement. The bond-length difference, Ar, between 
the end-members is 0.173 A and is ^ 0.14 A in the al- 
loy If we add the discrete displacements on the 
arsenic and metal sites we get Ar = 0.16(2) A. The 
bond-length difference can be obtained directly by fitting 
Gaussians to the first PDF peak in a model independent 
way What this modeling shows is that within the 
structure, most of the relaxation of local bonds occurs 
by arsenic moving off its site, but displacements of the 
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FIG. 7. Schematic representation of the discrete distortions 
to the zinc-blend structure obtained from the fitting. Frag- 
ments from the ideal structure (top panel) can be compared 
with the distorted (lower panel) structure. The view is a pro- 
jection along the [010] direction. The large circles are the As 
sites and the small circles the metal sites. The relative sizes of 
the circles reflect schematically the relative sizes of the con- 
tinuous-displacements on each site in the undoped material 
(top panel) and the alloy (bottom); however, their size has 
been exaggerated. 



metal atoms are also important. 

Ion channeling results give a precise determination of 



the mean-square displacement amplitude, 



(includ- 



ing static and dynamic components) perpendicular to the 
channeling direction. The results of Haga et al. give 
a (u^) for Gao.47Ino.53As of 0.017 at room tempera- 
ture. Based on the theoretical thermal amplitude of the 
end-member compounds being {u ) = 0.0121 A2, (and 
in good agreement with the value measured using the 
PDF j20j) they determined that the mean-square static 
displacements perpendicular to [100] were of magnitude 
(m^) — 0.005 A^. They saw a similar value in directions 
perpendicular to (110). This corresponds to a static root- 
mean-square displacement amplitude of 0.07 A. This is 
much smaller than the discrete displacements of 0.133 A 
that we observe. However, if we find the average dis- 
placement by adding the discrete displacements on the 
As and metal sublattices in quadrature and dividing by 
2 we get 0.068 A^ in good agreement with the ion chan- 
neling. Their work did not report which sublattice con- 
tributed most of the disorder; however, there is a sug- 



gestion from electron diffraction [Q, in agreement with 
theory [^|J^[l0|] , that the As sublattice is more disordered 
as we show directly from our measurement. 

Finally, we note that the actual displacement pattern 
on the arsenic site is expected to have (111) type dis- 
placements (as in our model) but also significant (100) 
type displacements [^^ 10 1. In fact, recent calculations 
indicate that the (100) displacements should be signif- 
icantly more pronounced than (111) displacements, es- 
pecially at room temperature []l0[. We are undertak- 
ing a more sophisticated modeling approach to explore 
this prediction. We tried a simple (100) displaced model, 
analogous to the one described here, but found it to ex- 
plain the data less successfully than the (111) model we 
described. We feel that this is a deficiency of the simple 
single-displacement-direction modeling rather than signi- 
fying that the displacement directions are actually (111) 
in the real alloy. The likely reason is that larger displace- 
ment amplitudes are required along (100) directions to 
satisfy the bond-length difference seen in the first PDF 
peak (actually times larger). With these large dis- 
crete displacements it is harder for the model to account 
for additional disorder in the data using enlarged ther- 
mal factors. A better fit in these imperfect models is ob- 
tained with smaller discrete displacements coupled with 
larger continuous-displacements. However, an improved 
fit should result from a more sophisticated model which 
includes both (100) and (111) type displacements. 



V. CONCLUSIONS 

From high real space resolution total and In differential 
PDFs of Ino.5Gao.5As alloy we conclude the following: 
In good agreement with earlier XAFS results [^ the Ga- 
As and In-As bonds do not take some compositionally 
averaged length but remain close to their natural lengths 
in Ino.5Gao.5As alloy. This bond- length mismatch brings 
about a considerable local disorder seen as a significant 
broadening of the next- nearest atomic-pair distributions. 
The positions and widths of the low and high-r peaks in 
both the total- and indium differential-PDFs are very 
well reproduced using a relaxed supercell model based 
on the Kirkwood potential with parameters taken from 
fits to the end-members of the alloy series. This suggests 
that this is a reasonable approach for generating the local 
structure of these alloys. 

A co-refinement of both the high-resolution total PDF 
and the chemically resolved indium differential PDF us- 
ing a simplified structural model was carried out. The 
arsenic sublattice contains most of the disorder in the 
structure as evidenced from both discrete atomic dis- 
placements in the model and enlarged thermal param- 
eters. However, small, but significant, displacements are 
evident on the metal sites and these have been quantified. 
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